This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.
Issue Howard Bayes
Time series 1999 - present 1977 - present
Bias in SWHS Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests. Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data. Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas. Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with. By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.
**Figure 1.**- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.


**Figure 2.**- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).


Results

**Figure X.**- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

**Figure 3.**- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 3.**- Total rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

**Figure 4.**- Black rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


And black rockfish releases…

**Figure 5.**- Black rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 6.**- Yellow rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 7.**- Yellow rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 8.**- DSR rockfish (excluding yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 9.**- DSR rockfish releases (including yelloweye) 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 11.**- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 12.**- Slope rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

**Figure 13.**- *Other* rockfish (non-pelagic and non-yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 14.**- *Other* rockfish (non-pelagic and non-yelloweye) releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



Total Biomass Removal Estimates

**Figure 15.**- Black rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.



**Figure 16.**- Yellow rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

**Figure 17.**- Pelagic rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 18.**- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 19.**- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


Model fit

Logbook residuals

**Figure 20.**- Residuals from logbook harvests.

Figure 20.- Residuals from logbook harvests.


SWHS residuals

**Figure 21.**- Residuals from SWHS harvests.

Figure 21.- Residuals from SWHS harvests.



**Figure 22.**- Residual of SWHS releases.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

**Figure 23.**- Mean percent of harvest by charter anglers.

Figure 23.- Mean percent of harvest by charter anglers.


When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

**Figure 24.**- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.


**Figure 25.**- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 26.**- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 27.**- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.


**Figure 28.**- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.


**Figure 29.**- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

**Figure 30.**- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.


Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

**Figure 31.**- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

**Figure 32.**- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

**Figure 33.**- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

**Figure 34.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

**Figure 35.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

**Figure 36.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.



P(slope|non-pelagic & non-yellowye) For release estimates

## NULL



Weight Fits

**Figure 37.**- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 38.**- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 39.**- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 40.**- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 41.**- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.5745856

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.11 while 0.89 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable

area_n

year

user

mean

median

lower

upper

Rhat

sd

skew

area

p_pelagic

11

1

1

0.455

0.453

0.220

0.709

1.011

0.128

0.016

CSEO

p_pelagic

11

1

2

0.391

0.384

0.173

0.647

1.015

0.124

0.056

CSEO

p_pelagic

11

10

2

0.389

0.380

0.179

0.645

1.011

0.122

0.074

CSEO

p_pelagic

11

11

1

0.449

0.449

0.209

0.707

1.010

0.131

0.000

CSEO

p_pelagic

11

11

2

0.391

0.383

0.179

0.635

1.012

0.122

0.065

CSEO

p_pelagic

11

12

2

0.390

0.380

0.178

0.643

1.012

0.121

0.083

CSEO

p_pelagic

11

13

1

0.457

0.454

0.225

0.714

1.018

0.129

0.023

CSEO

p_pelagic

11

13

2

0.388

0.379

0.174

0.641

1.011

0.123

0.073

CSEO

p_pelagic

11

14

1

0.465

0.460

0.240

0.695

1.019

0.121

0.041

CSEO

p_pelagic

11

16

2

0.396

0.392

0.174

0.646

1.011

0.123

0.033

CSEO

p_pelagic

11

17

1

0.459

0.457

0.232

0.696

1.012

0.121

0.017

CSEO

p_pelagic

11

18

1

0.464

0.464

0.239

0.701

1.010

0.122

0.000

CSEO

p_pelagic

11

18

2

0.397

0.391

0.183

0.650

1.011

0.121

0.049

CSEO

p_pelagic

11

2

1

0.452

0.447

0.213

0.709

1.011

0.129

0.039

CSEO

p_pelagic

11

2

2

0.392

0.388

0.172

0.640

1.019

0.124

0.032

CSEO

p_pelagic

11

21

1

0.472

0.467

0.242

0.715

1.011

0.122

0.041

CSEO

p_pelagic

11

3

1

0.455

0.451

0.222

0.704

1.013

0.128

0.031

CSEO

p_pelagic

11

3

2

0.392

0.385

0.177

0.649

1.011

0.124

0.057

CSEO

p_pelagic

11

4

1

0.454

0.452

0.221

0.694

1.011

0.124

0.016

CSEO

p_pelagic

11

4

2

0.389

0.385

0.178

0.639

1.013

0.122

0.033

CSEO

p_pelagic

11

5

1

0.457

0.456

0.220

0.712

1.014

0.128

0.008

CSEO

p_pelagic

11

5

2

0.390

0.383

0.172

0.638

1.016

0.124

0.057

CSEO

p_pelagic

11

6

1

0.456

0.452

0.222

0.702

1.010

0.126

0.032

CSEO

p_pelagic

11

7

1

0.454

0.449

0.221

0.705

1.017

0.128

0.039

CSEO

p_pelagic

11

7

2

0.391

0.385

0.173

0.641

1.013

0.120

0.050

CSEO

p_pelagic

11

8

1

0.456

0.453

0.220

0.703

1.011

0.128

0.024

CSEO

p_pelagic

11

9

1

0.454

0.452

0.213

0.712

1.021

0.127

0.016

CSEO

p_pelagic

11

9

2

0.389

0.383

0.172

0.644

1.014

0.123

0.049

CSEO

p_pelagic

13

1

1

0.236

0.212

0.052

0.549

1.241

0.133

0.180

NSEI

p_pelagic

13

1

2

0.093

0.076

0.016

0.269

1.235

0.066

0.259

NSEI

p_pelagic

13

10

1

0.268

0.246

0.078

0.559

1.176

0.127

0.173

NSEI

p_pelagic

13

10

2

0.108

0.092

0.028

0.280

1.146

0.066

0.241

NSEI

p_pelagic

13

11

1

0.276

0.257

0.085

0.569

1.126

0.127

0.149

NSEI

p_pelagic

13

11

2

0.110

0.094

0.026

0.291

1.119

0.068

0.235

NSEI

p_pelagic

13

12

1

0.280

0.263

0.087

0.563

1.113

0.124

0.137

NSEI

p_pelagic

13

12

2

0.113

0.099

0.029

0.276

1.119

0.064

0.220

NSEI

p_pelagic

13

13

1

0.286

0.267

0.093

0.575

1.116

0.125

0.152

NSEI

p_pelagic

13

13

2

0.117

0.103

0.030

0.287

1.114

0.067

0.209

NSEI

p_pelagic

13

14

1

0.277

0.262

0.096

0.542

1.100

0.116

0.129

NSEI

p_pelagic

13

14

2

0.118

0.103

0.034

0.277

1.096

0.063

0.236

NSEI

p_pelagic

13

15

1

0.257

0.239

0.089

0.527

1.080

0.112

0.161

NSEI

p_pelagic

13

15

2

0.115

0.103

0.035

0.260

1.073

0.059

0.204

NSEI

p_pelagic

13

16

1

0.294

0.279

0.108

0.551

1.060

0.115

0.130

NSEI

p_pelagic

13

16

2

0.127

0.114

0.038

0.295

1.058

0.067

0.194

NSEI

p_pelagic

13

17

1

0.313

0.299

0.126

0.575

1.063

0.116

0.120

NSEI

p_pelagic

13

17

2

0.135

0.121

0.044

0.295

1.048

0.066

0.211

NSEI

p_pelagic

13

18

1

0.324

0.315

0.133

0.574

1.046

0.114

0.079

NSEI

p_pelagic

13

18

2

0.140

0.126

0.044

0.313

1.047

0.070

0.200

NSEI

p_pelagic

13

19

1

0.324

0.312

0.134

0.569

1.033

0.114

0.105

NSEI

p_pelagic

13

19

2

0.146

0.133

0.049

0.324

1.035

0.071

0.184

NSEI

p_pelagic

13

2

1

0.242

0.218

0.056

0.548

1.264

0.131

0.183

NSEI

p_pelagic

13

2

2

0.094

0.077

0.018

0.269

1.227

0.067

0.253

NSEI

p_pelagic

13

20

1

0.343

0.331

0.151

0.586

1.024

0.115

0.105

NSEI

p_pelagic

13

20

2

0.155

0.142

0.051

0.325

1.027

0.071

0.182

NSEI

p_pelagic

13

21

1

0.363

0.352

0.166

0.615

1.023

0.116

0.095

NSEI

p_pelagic

13

21

2

0.165

0.152

0.061

0.335

1.014

0.073

0.179

NSEI

p_pelagic

13

22

2

0.168

0.157

0.057

0.337

1.013

0.074

0.149

NSEI

p_pelagic

13

3

1

0.243

0.219

0.059

0.538

1.226

0.129

0.187

NSEI

p_pelagic

13

3

2

0.096

0.079

0.019

0.274

1.231

0.067

0.255

NSEI

p_pelagic

13

4

1

0.246

0.223

0.061

0.547

1.239

0.132

0.175

NSEI

p_pelagic

13

4

2

0.097

0.080

0.019

0.263

1.212

0.066

0.258

NSEI

p_pelagic

13

5

1

0.247

0.224

0.063

0.552

1.232

0.128

0.180

NSEI

p_pelagic

13

5

2

0.098

0.082

0.020

0.264

1.187

0.064

0.248

NSEI

p_pelagic

13

6

1

0.250

0.228

0.066

0.546

1.212

0.129

0.170

NSEI

p_pelagic

13

6

2

0.101

0.086

0.022

0.268

1.186

0.067

0.224

NSEI

p_pelagic

13

7

1

0.254

0.232

0.069

0.554

1.167

0.128

0.172

NSEI

p_pelagic

13

7

2

0.102

0.086

0.021

0.271

1.186

0.068

0.236

NSEI

p_pelagic

13

8

1

0.261

0.242

0.071

0.559

1.194

0.129

0.147

NSEI

p_pelagic

13

8

2

0.102

0.086

0.023

0.277

1.172

0.067

0.239

NSEI

p_pelagic

13

9

1

0.265

0.245

0.078

0.551

1.176

0.125

0.160

NSEI

p_pelagic

13

9

2

0.105

0.088

0.023

0.277

1.144

0.067

0.255

NSEI

p_pelagic

14

21

1

0.438

0.433

0.214

0.691

1.012

0.123

0.041

NSEO

p_pelagic

14

6

1

0.428

0.422

0.204

0.687

1.010

0.123

0.049

NSEO

p_pelagic

3

12

2

0.191

0.178

0.067

0.395

1.010

0.085

0.154

PWSI

p_pelagic

3

2

2

0.189

0.176

0.065

0.391

1.010

0.085

0.152

PWSI

p_pelagic

4

7

2

0.557

0.559

0.301

0.794

1.011

0.129

-0.016

PWSO

p_black

8

37

1

0.381

0.379

0.317

0.449

1.065

0.042

0.047

afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable

area_n

year/species

mean

median

lower

upper

Rhat

sd

skew

year

area

beta2_pH

10

3

0.329

0.235

0.103

0.975

1.110

0.474

0.198

northeast

beta1_pH

11

1

3.252

3.171

2.489

4.387

1.138

0.487

0.166

CSEO

beta0_pelagic

11

-0.186

-0.167

-0.689

0.171

1.131

0.207

-0.092

CSEO

beta1_pelagic

11

1.231

1.194

0.711

2.118

1.187

0.323

0.114

CSEO

beta0_pH

13

4

0.735

0.816

-0.982

1.866

1.117

0.672

-0.121

NSEI

beta0_pelagic

13

-1.565

-1.465

-3.341

-0.232

1.845

0.830

-0.121

NSEI

beta1_pelagic

13

3.886

3.478

1.735

7.318

1.606

1.511

0.270

NSEI

beta2_pelagic

13

0.141

0.126

0.057

0.322

1.324

0.074

0.204

NSEI

beta1_pelagic

14

1.846

1.725

1.131

3.036

1.131

0.506

0.239

NSEO

beta2_yellow

14

-6.274

-4.748

-18.665

-0.376

1.203

4.981

-0.306

NSEO

beta2_yellow

15

-6.565

-5.503

-19.371

-0.277

1.134

5.093

-0.209

SSEI

beta2_yellow

2

-9.188

-7.428

-24.626

-0.602

1.130

6.774

-0.260

NG

beta0_pH

3

1

1.365

1.406

0.642

1.747

1.121

0.262

-0.157

PWSI

beta2_pH

8

1

0.442

0.278

0.121

2.000

1.205

0.687

0.239

afognak

beta2_pH

9

3

0.197

0.167

0.082

0.429

1.185

0.190

0.158

eastside

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.01

459.00

190.00

918.00

324.00

Prop. of total

0.71

0.29

0.74

0.26

Prop. median less than 1

0.28

0.24

0.32

0.44

Prop. median less than 10

0.36

0.60

0.52

0.77

Prop. median less than 50

0.48

0.80

0.65

0.85

Median of the median values

101.54

6.10

7.67

1.53

Mean of the median values

663.22

176.80

227.67

152.22

Mean skew of posteriors

0.24

0.29

0.27

0.25

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.1

40.00

6.00

51.00

32.00

Prop. of total

0.87

0.13

0.61

0.39

Prop. median less than 1

0.03

0.67

0.76

0.50

Prop. median less than 10

0.03

1.00

0.78

0.94

Prop. median less than 50

0.03

1.00

0.80

0.94

Median of the median values

745.76

0.38

0.01

1.07

Mean of the median values

656.44

2.68

61.55

6.00

Mean skew of posteriors

0.24

0.18

0.17

0.19

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable

Freq

beta0_pelagic

2

beta0_pH

2

beta1_pelagic

3

beta1_pH

1

beta2_pelagic

1

beta2_pH

3

beta2_yellow

3

Bp_ay

10

Bp_ayg

6

Bp_ayu

8

H_ayg

1

Hb_ay

10

Hb_ayg

9

Hb_ayu

6

Ho_ay

3

Ho_ayg

1

Ho_ayu

3

Hp_ay

11

Hp_ayg

11

Hp_ayu

8

Hs_ayu

4

Hy_ayu

1

p_pelagic

26

R_ayg

1

R_ayu

4

Rb_ayg

3

Rd_ayg

1

Rd_ayu

2

Ro_ay

9

Ro_ayg

3

Ro_ayu

25

Rp_ayg

1

Rp_ayu

4

Rs_ay

1

Rs_ayg

3

Rs_ayu

6

Ry_ay

7

Ry_ayg

3

Ry_ayu

20

Parameter estimates:

Summary Table of Parameter Estimates
Parameter mean sd Lower_CI Median Upper_CI
mu_bc_H[1] -0.015 0.070 -0.143 -0.018 0.133
mu_bc_H[2] 0.076 0.055 -0.024 0.073 0.192
mu_bc_H[3] -0.269 0.066 -0.400 -0.269 -0.137
mu_bc_H[4] -0.925 0.185 -1.312 -0.917 -0.577
mu_bc_H[5] 1.313 1.055 0.102 1.065 3.987
mu_bc_H[6] -2.407 0.320 -3.014 -2.413 -1.758
mu_bc_H[7] -0.455 0.100 -0.666 -0.453 -0.265
mu_bc_H[8] 0.266 0.344 -0.325 0.234 1.023
mu_bc_H[9] -0.244 0.126 -0.486 -0.244 0.013
mu_bc_H[10] 0.058 0.064 -0.061 0.056 0.190
mu_bc_H[11] -0.126 0.037 -0.202 -0.127 -0.055
mu_bc_H[12] -0.249 0.107 -0.482 -0.244 -0.044
mu_bc_H[13] -0.159 0.071 -0.296 -0.162 -0.018
mu_bc_H[14] -0.308 0.089 -0.490 -0.306 -0.141
mu_bc_H[15] -0.335 0.048 -0.426 -0.335 -0.242
mu_bc_H[16] -0.146 0.359 -0.762 -0.179 0.678
mu_bc_R[1] 1.042 0.120 0.808 1.043 1.283
mu_bc_R[2] 1.014 0.075 0.856 1.017 1.153
mu_bc_R[3] 0.223 0.089 0.042 0.224 0.391
mu_bc_R[4] 0.665 0.194 0.265 0.674 1.018
mu_bc_R[5] 1.674 0.369 0.960 1.680 2.409
mu_bc_R[6] -1.436 0.343 -2.084 -1.441 -0.748
mu_bc_R[7] 0.520 0.169 0.165 0.528 0.839
mu_bc_R[8] 0.647 0.148 0.357 0.646 0.942
mu_bc_R[9] 0.385 0.168 0.035 0.392 0.709
mu_bc_R[10] 0.895 0.103 0.697 0.892 1.098
mu_bc_R[11] 1.101 0.075 0.956 1.100 1.248
mu_bc_R[12] 0.873 0.162 0.550 0.875 1.198
mu_bc_R[13] 1.029 0.088 0.844 1.032 1.196
mu_bc_R[14] 0.859 0.123 0.618 0.859 1.100
mu_bc_R[15] 0.780 0.098 0.583 0.777 0.974
mu_bc_R[16] 1.270 0.115 1.035 1.273 1.484
tau_pH[1] 2.464 0.180 2.119 2.462 2.833
tau_pH[2] 1.566 0.144 1.301 1.559 1.867
tau_pH[3] 2.268 0.301 1.735 2.251 2.907
tau_pH[4] 1.124 0.175 0.818 1.109 1.500
tau_pH[5] 0.788 0.131 0.555 0.779 1.062
beta0_pH[1,1] 0.556 0.176 0.178 0.564 0.877
beta0_pH[2,1] 1.418 0.176 1.029 1.429 1.728
beta0_pH[3,1] 1.365 0.262 0.642 1.406 1.747
beta0_pH[4,1] 1.528 0.259 0.927 1.557 1.959
beta0_pH[5,1] -0.811 0.353 -1.572 -0.790 -0.192
beta0_pH[6,1] -0.544 0.444 -1.711 -0.438 0.046
beta0_pH[7,1] 0.752 0.011 0.732 0.752 0.773
beta0_pH[8,1] -0.694 0.330 -1.436 -0.665 -0.118
beta0_pH[9,1] -0.350 0.317 -1.144 -0.316 0.163
beta0_pH[10,1] 0.263 0.194 -0.192 0.269 0.631
beta0_pH[11,1] 0.107 0.306 -0.546 0.149 0.592
beta0_pH[12,1] 0.615 0.192 0.230 0.618 0.983
beta0_pH[13,1] -0.656 0.147 -0.936 -0.655 -0.371
beta0_pH[14,1] 0.012 0.172 -0.340 0.020 0.341
beta0_pH[15,1] 0.064 0.180 -0.331 0.078 0.376
beta0_pH[16,1] -0.874 0.189 -1.258 -0.869 -0.524
beta0_pH[1,2] 2.566 0.083 2.406 2.565 2.733
beta0_pH[2,2] 3.078 0.022 3.036 3.078 3.121
beta0_pH[3,2] 3.625 0.041 3.547 3.624 3.704
beta0_pH[4,2] 2.936 0.019 2.900 2.936 2.973
beta0_pH[5,2] 4.568 1.461 2.788 4.269 8.330
beta0_pH[6,2] 3.327 0.161 3.035 3.321 3.656
beta0_pH[7,2] 2.026 0.050 1.931 2.025 2.122
beta0_pH[8,2] 2.824 0.065 2.703 2.823 2.958
beta0_pH[9,2] 3.963 0.210 3.561 3.958 4.387
beta0_pH[10,2] 3.998 0.129 3.752 3.993 4.262
beta0_pH[11,2] 1.889 0.009 1.871 1.889 1.906
beta0_pH[12,2] 1.804 0.070 1.673 1.802 1.948
beta0_pH[13,2] 2.496 0.026 2.446 2.496 2.547
beta0_pH[14,2] 1.806 0.024 1.760 1.807 1.854
beta0_pH[15,2] 2.550 0.017 2.516 2.549 2.582
beta0_pH[16,2] 2.776 0.024 2.730 2.776 2.822
beta0_pH[1,3] 0.718 0.566 -0.771 0.872 1.374
beta0_pH[2,3] 2.126 0.018 2.092 2.126 2.161
beta0_pH[3,3] 2.689 0.024 2.642 2.689 2.735
beta0_pH[4,3] 2.981 0.034 2.914 2.981 3.049
beta0_pH[5,3] 3.722 2.649 0.928 3.020 11.045
beta0_pH[6,3] 1.245 0.321 0.675 1.228 1.965
beta0_pH[7,3] 0.855 0.050 0.760 0.854 0.952
beta0_pH[8,3] 0.283 0.125 0.018 0.290 0.505
beta0_pH[9,3] -2.742 1.068 -5.193 -2.590 -0.979
beta0_pH[10,3] -0.827 0.844 -2.741 -0.743 0.593
beta0_pH[11,3] 0.037 11.484 -22.266 0.016 22.574
beta0_pH[12,3] 0.053 11.361 -22.389 0.130 22.061
beta0_pH[13,3] 0.108 11.507 -22.806 0.257 22.574
beta0_pH[14,3] 0.001 11.624 -23.068 0.095 22.375
beta0_pH[15,3] 0.099 11.270 -21.706 0.200 21.599
beta0_pH[16,3] 0.135 11.362 -22.357 0.015 21.805
beta0_pH[11,4] 0.931 0.763 -1.331 1.023 2.115
beta0_pH[12,4] 0.820 0.816 -1.041 0.843 2.296
beta0_pH[13,4] 0.735 0.672 -0.982 0.816 1.866
beta0_pH[14,4] 0.291 0.578 -1.091 0.360 1.224
beta0_pH[15,4] 0.577 0.618 -1.048 0.646 1.566
beta0_pH[16,4] 0.981 0.739 -0.697 1.004 2.323
beta0_pH[11,5] 1.671 0.088 1.513 1.665 1.857
beta0_pH[12,5] -1.410 0.450 -2.334 -1.386 -0.569
beta0_pH[13,5] 2.233 0.259 1.811 2.205 2.817
beta0_pH[14,5] 1.894 0.304 1.379 1.862 2.561
beta0_pH[15,5] 1.109 0.065 0.996 1.105 1.243
beta0_pH[16,5] -0.737 0.051 -0.836 -0.737 -0.636
beta1_pH[1,1] 2.575 0.480 1.824 2.514 3.813
beta1_pH[2,1] 1.606 0.344 1.019 1.569 2.360
beta1_pH[3,1] 1.858 0.480 1.139 1.790 3.221
beta1_pH[4,1] 2.028 0.535 1.255 1.914 3.565
beta1_pH[5,1] 2.338 0.478 1.560 2.286 3.410
beta1_pH[6,1] 2.264 0.903 1.089 2.020 4.551
beta1_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,1] 2.701 0.882 1.369 2.524 4.566
beta1_pH[9,1] 2.205 0.730 1.264 2.024 4.267
beta1_pH[10,1] 2.190 0.404 1.528 2.143 3.250
beta1_pH[11,1] 3.252 0.487 2.489 3.171 4.387
beta1_pH[12,1] 2.921 0.287 2.356 2.923 3.482
beta1_pH[13,1] 3.455 0.244 2.992 3.449 3.935
beta1_pH[14,1] 2.846 0.319 2.229 2.838 3.485
beta1_pH[15,1] 2.385 0.261 1.920 2.365 2.944
beta1_pH[16,1] 4.037 0.271 3.526 4.032 4.582
beta1_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[1,3] 3.403 1.688 1.224 2.926 7.396
beta1_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,3] 4.888 2.303 2.249 4.108 11.268
beta1_pH[9,3] 8.158 2.544 4.124 7.721 13.859
beta1_pH[10,3] 7.012 2.601 3.113 6.539 12.801
beta1_pH[11,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,4] 2.240 0.786 0.880 2.172 4.411
beta1_pH[12,4] 3.339 0.823 1.839 3.298 5.262
beta1_pH[13,4] 1.741 0.709 0.433 1.698 3.489
beta1_pH[14,4] 2.364 0.660 1.251 2.314 3.884
beta1_pH[15,4] 1.865 0.671 0.731 1.815 3.557
beta1_pH[16,4] 3.880 0.790 2.425 3.859 5.649
beta1_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,1] 1.228 0.825 0.338 0.937 3.303
beta2_pH[2,1] 1.541 0.952 0.295 1.314 3.408
beta2_pH[3,1] 1.304 0.833 0.182 1.148 2.899
beta2_pH[4,1] 1.014 0.757 0.175 0.750 2.852
beta2_pH[5,1] 2.039 1.559 0.259 1.684 6.148
beta2_pH[6,1] 0.863 1.260 0.104 0.377 4.261
beta2_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,1] 0.442 0.687 0.121 0.278 2.000
beta2_pH[9,1] 0.836 1.177 0.132 0.472 4.339
beta2_pH[10,1] 1.141 1.124 0.240 0.801 4.454
beta2_pH[11,1] 0.466 0.261 0.200 0.399 1.075
beta2_pH[12,1] 1.526 0.798 0.609 1.265 3.660
beta2_pH[13,1] 1.023 0.411 0.550 0.928 2.171
beta2_pH[14,1] 1.097 0.629 0.450 0.913 3.103
beta2_pH[15,1] 1.028 0.595 0.399 0.874 2.891
beta2_pH[16,1] 1.052 0.345 0.624 0.990 1.909
beta2_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,3] 2.518 4.348 0.115 0.524 15.912
beta2_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,3] 2.436 2.014 0.564 1.597 8.183
beta2_pH[9,3] 0.197 0.190 0.082 0.167 0.429
beta2_pH[10,3] 0.329 0.474 0.103 0.235 0.975
beta2_pH[11,3] -8.397 5.982 -22.684 -7.170 -0.424
beta2_pH[12,3] -8.429 5.993 -22.654 -7.215 -0.396
beta2_pH[13,3] -8.371 6.015 -22.575 -7.136 -0.382
beta2_pH[14,3] -8.359 5.952 -22.633 -7.151 -0.400
beta2_pH[15,3] -8.349 5.929 -22.463 -7.300 -0.404
beta2_pH[16,3] -8.370 5.941 -22.789 -7.216 -0.478
beta2_pH[11,4] -5.188 4.274 -16.623 -3.892 -0.641
beta2_pH[12,4] -5.767 4.199 -16.782 -4.517 -1.011
beta2_pH[13,4] -5.337 4.255 -16.622 -3.898 -0.847
beta2_pH[14,4] -5.370 4.220 -16.344 -3.868 -0.969
beta2_pH[15,4] -5.439 4.154 -16.773 -4.148 -0.826
beta2_pH[16,4] -5.530 3.971 -16.322 -4.212 -1.295
beta2_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta3_pH[1,1] 34.990 1.166 32.915 34.969 37.655
beta3_pH[2,1] 32.728 1.277 30.357 32.698 35.272
beta3_pH[3,1] 32.986 1.627 29.500 33.030 36.217
beta3_pH[4,1] 32.841 1.960 29.620 32.653 37.478
beta3_pH[5,1] 27.637 1.144 26.211 27.471 30.966
beta3_pH[6,1] 33.546 3.342 27.194 33.103 41.303
beta3_pH[7,1] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,1] 35.857 3.696 29.043 36.123 42.690
beta3_pH[9,1] 32.732 3.084 27.939 32.156 40.341
beta3_pH[10,1] 32.658 1.359 30.377 32.536 35.515
beta3_pH[11,1] 31.229 1.016 29.307 31.239 33.207
beta3_pH[12,1] 30.776 0.469 29.807 30.797 31.615
beta3_pH[13,1] 31.413 0.538 30.400 31.399 32.510
beta3_pH[14,1] 32.497 0.739 31.243 32.450 34.031
beta3_pH[15,1] 30.922 0.726 29.469 30.949 32.301
beta3_pH[16,1] 30.114 0.355 29.386 30.126 30.797
beta3_pH[1,2] 34.411 5.471 25.551 34.267 43.442
beta3_pH[2,2] 34.746 5.464 25.468 35.019 43.448
beta3_pH[3,2] 30.502 7.257 18.579 30.458 42.435
beta3_pH[4,2] 30.451 7.171 18.591 30.280 42.384
beta3_pH[5,2] 30.321 7.228 18.534 30.275 42.355
beta3_pH[6,2] 30.398 7.330 18.654 30.288 42.407
beta3_pH[7,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[9,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[10,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[11,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[12,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[13,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[14,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[15,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[16,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[1,3] 39.817 2.527 33.627 40.176 43.708
beta3_pH[2,3] 34.512 5.522 25.442 34.722 43.500
beta3_pH[3,3] 30.484 7.269 18.637 30.539 42.381
beta3_pH[4,3] 30.287 7.229 18.571 30.262 42.331
beta3_pH[5,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[6,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[7,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,3] 42.134 0.831 40.912 41.872 43.841
beta3_pH[9,3] 35.876 1.622 34.521 35.287 40.712
beta3_pH[10,3] 37.690 2.440 34.634 37.071 43.239
beta3_pH[11,3] 38.498 2.589 34.239 38.513 42.751
beta3_pH[12,3] 38.506 2.606 34.207 38.512 42.740
beta3_pH[13,3] 38.542 2.584 34.265 38.574 42.731
beta3_pH[14,3] 38.495 2.595 34.215 38.579 42.754
beta3_pH[15,3] 38.437 2.575 34.215 38.374 42.797
beta3_pH[16,3] 38.448 2.594 34.265 38.446 42.791
beta3_pH[11,4] 40.811 1.181 38.877 40.700 43.249
beta3_pH[12,4] 41.923 0.789 40.262 42.185 42.936
beta3_pH[13,4] 41.402 1.074 39.280 41.364 43.598
beta3_pH[14,4] 40.399 0.744 38.912 40.428 41.810
beta3_pH[15,4] 41.197 0.890 39.597 41.121 43.155
beta3_pH[16,4] 41.753 0.426 40.915 41.813 42.566
beta3_pH[11,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[12,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[13,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[14,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[15,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[16,5] 25.000 0.000 25.000 25.000 25.000
beta4_pH[1,1] 0.262 0.021 0.220 0.262 0.304
beta4_pH[2,1] 0.211 0.005 0.201 0.211 0.221
beta4_pH[3,1] 0.172 0.007 0.158 0.172 0.187
beta4_pH[4,1] 0.294 0.011 0.272 0.294 0.315
beta4_pH[5,1] 0.434 0.335 0.089 0.390 0.981
beta4_pH[6,1] 0.784 0.394 0.326 0.735 1.459
beta4_pH[7,1] 0.853 0.343 0.342 0.801 1.693
beta4_pH[8,1] 0.155 0.070 0.042 0.148 0.314
beta4_pH[9,1] 0.217 0.086 0.063 0.211 0.392
beta4_pH[10,1] 0.001 0.001 0.000 0.001 0.003
beta4_pH[11,1] 0.445 0.012 0.423 0.445 0.469
beta4_pH[12,1] 0.641 0.023 0.597 0.642 0.684
beta4_pH[13,1] 0.630 0.014 0.604 0.630 0.657
beta4_pH[14,1] 0.669 0.031 0.609 0.668 0.730
beta4_pH[15,1] 0.798 0.021 0.758 0.797 0.839
beta4_pH[16,1] 1.214 0.041 1.137 1.213 1.293
beta4_pH[1,2] 0.298 0.098 0.122 0.293 0.503
beta4_pH[2,2] 0.514 0.022 0.470 0.513 0.558
beta4_pH[3,2] 0.446 0.012 0.423 0.446 0.470
beta4_pH[4,2] 0.655 0.038 0.583 0.655 0.731
beta4_pH[5,2] 0.711 0.723 0.109 0.533 2.367
beta4_pH[6,2] 0.907 0.756 0.203 0.743 2.606
beta4_pH[7,2] 0.896 0.508 0.339 0.811 1.998
beta4_pH[8,2] 0.848 0.202 0.496 0.835 1.292
beta4_pH[9,2] 1.163 1.107 0.305 0.884 3.741
beta4_pH[10,2] 0.458 0.052 0.361 0.455 0.565
beta4_pH[11,2] 1.268 0.050 1.171 1.267 1.368
beta4_pH[12,2] 0.605 0.110 0.400 0.603 0.843
beta4_pH[13,2] 0.949 0.045 0.860 0.948 1.035
beta4_pH[14,2] 0.935 0.063 0.813 0.933 1.060
beta4_pH[15,2] 1.010 0.026 0.960 1.010 1.061
beta4_pH[16,2] 1.388 0.076 1.240 1.386 1.538
beta4_pH[1,3] 0.241 0.076 0.108 0.236 0.407
beta4_pH[2,3] 0.050 0.017 0.020 0.049 0.083
beta4_pH[3,3] 0.052 0.010 0.033 0.052 0.072
beta4_pH[4,3] 0.167 0.031 0.108 0.167 0.229
beta4_pH[5,3] 0.760 0.823 0.127 0.560 2.721
beta4_pH[6,3] 0.879 0.864 0.194 0.709 2.542
beta4_pH[7,3] 0.902 0.462 0.347 0.813 2.014
beta4_pH[8,3] 0.388 0.189 0.097 0.362 0.844
beta4_pH[9,3] 0.922 0.926 0.176 0.715 2.845
beta4_pH[10,3] 0.054 0.038 0.003 0.048 0.146
beta4_pH[11,3] 0.618 0.257 0.262 0.572 1.277
beta4_pH[12,3] 0.707 0.232 0.377 0.672 1.277
beta4_pH[13,3] 0.722 0.193 0.416 0.697 1.183
beta4_pH[14,3] 0.860 0.217 0.510 0.834 1.353
beta4_pH[15,3] 0.823 0.243 0.426 0.795 1.360
beta4_pH[16,3] 0.842 0.389 0.317 0.764 1.847
beta4_pH[1,4] 0.631 0.362 0.207 0.555 1.559
beta4_pH[2,4] 0.820 1.050 0.116 0.571 3.018
beta4_pH[3,4] 0.880 1.016 0.114 0.599 3.335
beta4_pH[4,4] 0.735 0.496 0.199 0.618 1.941
beta4_pH[5,4] 0.777 1.030 0.120 0.569 2.598
beta4_pH[6,4] 0.921 1.038 0.201 0.728 2.795
beta4_pH[7,4] 0.907 0.477 0.360 0.813 1.996
beta4_pH[8,4] 0.730 0.687 0.138 0.576 2.176
beta4_pH[9,4] 0.793 0.922 0.122 0.561 2.907
beta4_pH[10,4] 1.390 2.860 0.008 0.384 10.257
beta4_pH[11,4] 0.396 0.016 0.365 0.396 0.427
beta4_pH[12,4] 0.551 0.053 0.450 0.549 0.661
beta4_pH[13,4] 0.836 0.040 0.760 0.835 0.916
beta4_pH[14,4] 0.863 0.071 0.731 0.859 1.010
beta4_pH[15,4] 0.626 0.017 0.593 0.625 0.661
beta4_pH[16,4] 0.647 0.039 0.576 0.646 0.727
beta4_pH[1,5] 0.639 0.345 0.198 0.567 1.447
beta4_pH[2,5] 0.820 0.885 0.112 0.586 2.968
beta4_pH[3,5] 0.860 0.929 0.105 0.583 3.399
beta4_pH[4,5] 0.726 0.501 0.208 0.613 1.840
beta4_pH[5,5] 0.762 0.890 0.117 0.552 2.679
beta4_pH[6,5] 0.896 0.739 0.202 0.726 2.709
beta4_pH[7,5] 0.902 0.529 0.343 0.816 1.870
beta4_pH[8,5] 0.722 0.589 0.152 0.586 2.201
beta4_pH[9,5] 0.785 0.976 0.113 0.558 2.710
beta4_pH[10,5] 1.475 2.953 0.008 0.391 11.316
beta4_pH[11,5] 0.361 0.075 0.218 0.360 0.512
beta4_pH[12,5] 0.694 0.205 0.364 0.667 1.180
beta4_pH[13,5] 0.973 0.111 0.769 0.966 1.210
beta4_pH[14,5] 0.842 0.159 0.561 0.831 1.184
beta4_pH[15,5] 1.311 0.120 1.091 1.303 1.564
beta4_pH[16,5] 0.922 0.379 0.376 0.860 1.845
beta0_pelagic[1] 2.268 0.036 2.202 2.267 2.339
beta0_pelagic[2] 1.538 0.007 1.524 1.538 1.552
beta0_pelagic[3] 0.180 0.222 -0.327 0.207 0.537
beta0_pelagic[4] 0.231 0.211 -0.189 0.233 0.640
beta0_pelagic[5] 1.445 0.160 1.144 1.443 1.762
beta0_pelagic[6] 1.356 0.061 1.238 1.354 1.479
beta0_pelagic[7] 1.470 0.038 1.398 1.468 1.547
beta0_pelagic[8] 1.672 0.044 1.592 1.671 1.760
beta0_pelagic[9] 2.775 0.064 2.655 2.774 2.899
beta0_pelagic[10] 2.603 0.031 2.543 2.602 2.666
beta0_pelagic[11] -0.186 0.207 -0.689 -0.167 0.171
beta0_pelagic[12] 1.801 0.035 1.735 1.800 1.871
beta0_pelagic[13] -1.565 0.830 -3.341 -1.465 -0.232
beta0_pelagic[14] -0.284 0.181 -0.697 -0.267 0.014
beta0_pelagic[15] -0.211 0.039 -0.327 -0.195 -0.179
beta0_pelagic[16] -0.445 0.350 -1.282 -0.350 0.011
beta1_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[3] 0.763 0.242 0.348 0.741 1.314
beta1_pelagic[4] 0.820 0.231 0.397 0.815 1.282
beta1_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[11] 1.231 0.323 0.711 1.194 2.118
beta1_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[13] 3.886 1.511 1.735 3.478 7.318
beta1_pelagic[14] 1.846 0.506 1.131 1.725 3.036
beta1_pelagic[15] 0.079 0.138 0.000 0.013 0.527
beta1_pelagic[16] 1.511 0.493 0.798 1.428 2.628
beta2_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[3] 6.230 5.971 0.209 4.487 21.044
beta2_pelagic[4] 7.649 6.040 0.459 6.341 22.243
beta2_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[11] 1.626 1.487 0.191 1.114 5.687
beta2_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[13] 0.141 0.074 0.057 0.126 0.322
beta2_pelagic[14] 0.840 0.842 0.188 0.576 3.366
beta2_pelagic[15] 2.588 1.915 0.151 2.187 7.231
beta2_pelagic[16] 0.739 0.920 0.140 0.390 3.907
beta3_pelagic[1] 28.600 3.893 22.298 28.479 35.098
beta3_pelagic[2] 28.765 3.882 22.323 28.797 35.125
beta3_pelagic[3] 27.577 2.704 23.248 27.509 32.257
beta3_pelagic[4] 26.268 1.545 23.580 26.298 30.428
beta3_pelagic[5] 29.368 6.665 18.525 29.324 40.422
beta3_pelagic[6] 29.585 6.603 18.620 29.365 40.482
beta3_pelagic[7] 29.204 6.629 18.458 29.075 40.232
beta3_pelagic[8] 29.405 6.595 18.574 29.372 40.387
beta3_pelagic[9] 29.576 6.656 18.551 29.726 40.456
beta3_pelagic[10] 29.407 6.558 18.631 29.392 40.458
beta3_pelagic[11] 30.154 1.249 27.864 30.053 32.989
beta3_pelagic[12] 30.000 0.000 30.000 30.000 30.000
beta3_pelagic[13] 30.972 4.950 23.501 30.390 41.067
beta3_pelagic[14] 34.452 1.940 31.486 34.147 39.318
beta3_pelagic[15] 33.276 5.805 23.607 33.593 41.703
beta3_pelagic[16] 30.982 2.326 26.114 31.082 35.695
beta0_yellow[1] -0.483 0.081 -0.641 -0.483 -0.326
beta0_yellow[2] 0.544 0.077 0.380 0.550 0.673
beta0_yellow[3] -0.212 0.129 -0.463 -0.217 0.013
beta0_yellow[4] 0.888 0.269 0.086 0.969 1.189
beta0_yellow[5] 0.236 0.323 -0.437 0.254 0.806
beta0_yellow[6] 0.819 0.107 0.611 0.820 1.025
beta0_yellow[7] 1.012 0.099 0.813 1.013 1.206
beta0_yellow[8] 1.020 0.093 0.836 1.021 1.193
beta0_yellow[9] 0.441 0.100 0.242 0.441 0.637
beta0_yellow[10] 0.537 0.070 0.400 0.536 0.677
beta0_yellow[11] -0.202 0.121 -0.455 -0.192 0.015
beta0_yellow[12] -1.083 0.093 -1.264 -1.082 -0.902
beta0_yellow[13] -1.335 0.177 -1.737 -1.315 -1.051
beta0_yellow[14] -0.243 0.144 -0.552 -0.230 0.001
beta0_yellow[15] -1.239 0.128 -1.532 -1.228 -1.035
beta0_yellow[16] -0.260 0.092 -0.444 -0.257 -0.086
beta1_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[2] 1.015 0.222 0.601 1.009 1.456
beta1_yellow[3] 0.467 0.239 0.020 0.479 0.909
beta1_yellow[4] 1.331 0.593 0.647 1.185 3.272
beta1_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[11] 0.294 0.369 0.001 0.106 1.268
beta1_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[13] 1.942 1.709 0.033 1.133 5.944
beta1_yellow[14] 0.714 0.967 0.012 0.435 3.768
beta1_yellow[15] 0.162 0.207 0.011 0.077 0.727
beta1_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[2] -9.188 6.774 -24.626 -7.428 -0.602
beta2_yellow[3] -7.645 5.827 -21.785 -6.305 -0.438
beta2_yellow[4] -5.649 5.893 -20.607 -3.690 -0.121
beta2_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[11] -7.937 5.917 -22.288 -6.790 -0.415
beta2_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[13] -5.454 5.361 -18.926 -3.472 -0.320
beta2_yellow[14] -6.274 4.981 -18.665 -4.748 -0.376
beta2_yellow[15] -6.565 5.093 -19.371 -5.503 -0.277
beta2_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta3_yellow[1] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[2] 28.840 0.995 27.454 28.667 31.617
beta3_yellow[3] 32.424 3.314 22.743 32.625 38.571
beta3_yellow[4] 29.953 3.361 24.425 28.680 35.991
beta3_yellow[5] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[6] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[7] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[8] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[9] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[10] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[11] 35.032 3.700 30.344 34.014 43.298
beta3_yellow[12] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[13] 32.394 2.798 29.092 32.047 38.557
beta3_yellow[14] 33.792 3.537 29.189 33.356 42.783
beta3_yellow[15] 36.405 3.869 30.259 36.247 43.501
beta3_yellow[16] 35.000 0.000 35.000 35.000 35.000
beta0_black[1] 0.013 0.077 -0.139 0.012 0.160
beta0_black[2] 1.855 0.032 1.795 1.855 1.916
beta0_black[3] 1.239 0.037 1.167 1.239 1.311
beta0_black[4] 2.425 0.047 2.335 2.426 2.518
beta0_black[5] 5.245 2.414 1.919 4.787 11.410
beta0_black[6] 5.230 2.155 2.593 4.653 10.999
beta0_black[7] 4.163 2.182 1.619 3.594 10.157
beta0_black[8] 0.884 0.217 0.488 0.880 1.338
beta0_black[9] 2.820 0.268 2.299 2.813 3.361
beta0_black[10] 1.288 0.051 1.188 1.288 1.389
beta0_black[11] 2.979 0.051 2.879 2.978 3.078
beta0_black[12] 4.917 0.196 4.530 4.915 5.310
beta0_black[13] 0.441 0.337 0.019 0.358 1.236
beta0_black[14] 2.853 0.106 2.634 2.855 3.067
beta0_black[15] 0.722 0.052 0.620 0.723 0.820
beta0_black[16] 3.772 0.075 3.627 3.769 3.921
beta2_black[1] 0.000 0.000 0.000 0.000 0.000
beta2_black[2] 0.000 0.000 0.000 0.000 0.000
beta2_black[3] 0.000 0.000 0.000 0.000 0.000
beta2_black[4] 0.000 0.000 0.000 0.000 0.000
beta2_black[5] 0.000 0.000 0.000 0.000 0.000
beta2_black[6] 0.000 0.000 0.000 0.000 0.000
beta2_black[7] 0.000 0.000 0.000 0.000 0.000
beta2_black[8] 0.000 0.000 0.000 0.000 0.000
beta2_black[9] 0.000 0.000 0.000 0.000 0.000
beta2_black[10] 0.000 0.000 0.000 0.000 0.000
beta2_black[11] 0.000 0.000 0.000 0.000 0.000
beta2_black[12] 0.000 0.000 0.000 0.000 0.000
beta2_black[13] -2.540 1.794 -6.885 -2.132 -0.254
beta2_black[14] 0.000 0.000 0.000 0.000 0.000
beta2_black[15] 0.000 0.000 0.000 0.000 0.000
beta2_black[16] 0.000 0.000 0.000 0.000 0.000
beta3_black[1] 25.000 0.000 25.000 25.000 25.000
beta3_black[2] 25.000 0.000 25.000 25.000 25.000
beta3_black[3] 25.000 0.000 25.000 25.000 25.000
beta3_black[4] 25.000 0.000 25.000 25.000 25.000
beta3_black[5] 25.000 0.000 25.000 25.000 25.000
beta3_black[6] 25.000 0.000 25.000 25.000 25.000
beta3_black[7] 25.000 0.000 25.000 25.000 25.000
beta3_black[8] 25.000 0.000 25.000 25.000 25.000
beta3_black[9] 25.000 0.000 25.000 25.000 25.000
beta3_black[10] 25.000 0.000 25.000 25.000 25.000
beta3_black[11] 25.000 0.000 25.000 25.000 25.000
beta3_black[12] 25.000 0.000 25.000 25.000 25.000
beta3_black[13] 37.297 3.437 26.139 38.108 40.807
beta3_black[14] 25.000 0.000 25.000 25.000 25.000
beta3_black[15] 25.000 0.000 25.000 25.000 25.000
beta3_black[16] 25.000 0.000 25.000 25.000 25.000
beta4_black[1] -0.325 0.092 -0.508 -0.322 -0.151
beta4_black[2] 0.258 0.052 0.157 0.259 0.362
beta4_black[3] -0.721 0.063 -0.847 -0.720 -0.595
beta4_black[4] 0.628 0.128 0.382 0.626 0.891
beta4_black[5] 0.735 1.484 -1.237 0.461 4.048
beta4_black[6] 0.865 1.507 -0.952 0.551 4.388
beta4_black[7] 0.777 1.341 -0.935 0.497 4.160
beta4_black[8] 0.071 0.288 -0.546 0.078 0.611
beta4_black[9] 1.040 0.819 -0.076 0.862 3.062
beta4_black[10] 0.074 0.067 -0.053 0.074 0.211
beta4_black[11] -0.610 0.064 -0.732 -0.610 -0.484
beta4_black[12] 1.239 0.284 0.718 1.227 1.810
beta4_black[13] -1.001 0.077 -1.155 -1.002 -0.854
beta4_black[14] 0.221 0.185 -0.143 0.220 0.584
beta4_black[15] 0.001 0.047 -0.090 0.001 0.092
beta4_black[16] -0.418 0.116 -0.641 -0.417 -0.187
beta0_dsr[11] -0.337 0.010 -0.356 -0.337 -0.319
beta0_dsr[12] 1.421 0.044 1.335 1.421 1.510
beta0_dsr[13] 0.130 0.022 0.087 0.131 0.172
beta0_dsr[14] -0.135 0.017 -0.169 -0.135 -0.102
beta0_dsr[15] 0.184 0.014 0.156 0.184 0.212
beta0_dsr[16] 0.208 0.015 0.178 0.208 0.237
beta1_dsr[11] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[12] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[13] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[14] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[15] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[16] 0.000 0.000 0.000 0.000 0.000
beta4_dsr[11] 0.548 0.035 0.482 0.548 0.618
beta4_dsr[12] 0.824 0.233 0.388 0.822 1.286
beta4_dsr[13] -0.048 0.041 -0.129 -0.047 0.033
beta4_dsr[14] -0.207 0.071 -0.347 -0.207 -0.065
beta4_dsr[15] 0.237 0.021 0.195 0.237 0.279
beta4_dsr[16] -0.099 0.045 -0.190 -0.099 -0.010
beta0_slope[11] -2.201 0.021 -2.241 -2.201 -2.160
beta0_slope[12] -5.439 0.225 -5.880 -5.437 -5.009
beta0_slope[13] -1.456 0.034 -1.523 -1.455 -1.390
beta0_slope[14] -2.692 0.050 -2.793 -2.691 -2.597
beta0_slope[15] -1.567 0.022 -1.610 -1.567 -1.524
beta0_slope[16] -3.015 0.045 -3.104 -3.015 -2.925
beta1_slope[11] 0.000 0.000 0.000 0.000 0.000
beta1_slope[12] 0.000 0.000 0.000 0.000 0.000
beta1_slope[13] 0.000 0.000 0.000 0.000 0.000
beta1_slope[14] 0.000 0.000 0.000 0.000 0.000
beta1_slope[15] 0.000 0.000 0.000 0.000 0.000
beta1_slope[16] 0.000 0.000 0.000 0.000 0.000
beta4_slope[11] -0.343 0.075 -0.494 -0.343 -0.197
beta4_slope[12] -4.643 0.819 -6.341 -4.590 -3.183
beta4_slope[13] 0.282 0.054 0.171 0.283 0.384
beta4_slope[14] -0.179 0.152 -0.479 -0.177 0.115
beta4_slope[15] -0.386 0.033 -0.449 -0.386 -0.322
beta4_slope[16] 0.221 0.104 0.016 0.221 0.430
beta5_slope[11] 3.076 0.109 2.873 3.073 3.301
beta5_slope[12] 1.992 0.391 1.260 1.975 2.799
beta5_slope[13] 0.931 0.081 0.766 0.933 1.087
beta5_slope[14] 3.780 0.368 3.101 3.756 4.554
beta5_slope[15] 1.387 0.072 1.248 1.385 1.529
beta5_slope[16] 2.684 0.235 2.243 2.671 3.179
tau_comp[1] 3.757 0.292 3.211 3.746 4.348
tau_comp[2] 2.852 0.251 2.372 2.847 3.360
tau_comp[3] 0.872 0.073 0.735 0.870 1.021
tau_comp[4] 1.274 0.194 0.933 1.258 1.702
tau_comp[5] 1.209 0.181 0.893 1.194 1.599
sigma_H[1] 0.209 0.048 0.122 0.207 0.309
sigma_H[2] 0.174 0.029 0.125 0.172 0.237
sigma_H[3] 0.188 0.040 0.117 0.185 0.276
sigma_H[4] 0.460 0.079 0.329 0.452 0.636
sigma_H[5] 0.967 0.201 0.594 0.954 1.396
sigma_H[6] 0.550 0.196 0.173 0.543 0.950
sigma_H[7] 0.310 0.064 0.212 0.302 0.456
sigma_H[8] 0.423 0.091 0.293 0.412 0.627
sigma_H[9] 0.513 0.118 0.331 0.494 0.773
sigma_H[10] 0.210 0.039 0.143 0.206 0.298
sigma_H[11] 0.266 0.044 0.192 0.263 0.366
sigma_H[12] 0.380 0.148 0.195 0.336 0.712
sigma_H[13] 0.218 0.036 0.155 0.215 0.297
sigma_H[14] 0.499 0.089 0.343 0.493 0.686
sigma_H[15] 0.235 0.038 0.171 0.232 0.321
sigma_H[16] 0.213 0.040 0.146 0.209 0.301
lambda_H[1] 2.331 3.631 0.084 1.075 11.572
lambda_H[2] 5.070 4.735 0.485 3.684 17.565
lambda_H[3] 5.437 8.082 0.277 2.703 27.979
lambda_H[4] 0.006 0.004 0.001 0.005 0.017
lambda_H[5] 2.635 6.306 0.028 0.678 19.775
lambda_H[6] 7.474 13.981 0.008 0.494 43.111
lambda_H[7] 0.012 0.008 0.002 0.010 0.033
lambda_H[8] 8.010 9.594 0.113 4.608 36.063
lambda_H[9] 0.015 0.010 0.003 0.013 0.039
lambda_H[10] 0.219 0.329 0.026 0.150 0.784
lambda_H[11] 0.233 0.326 0.011 0.105 1.144
lambda_H[12] 3.602 4.897 0.158 2.016 16.917
lambda_H[13] 3.192 3.107 0.255 2.285 11.324
lambda_H[14] 3.338 3.972 0.232 2.098 14.350
lambda_H[15] 0.019 0.020 0.003 0.014 0.061
lambda_H[16] 0.463 0.710 0.030 0.243 2.209
mu_lambda_H[1] 4.030 1.897 0.980 3.785 8.105
mu_lambda_H[2] 3.748 1.915 0.536 3.597 7.841
mu_lambda_H[3] 3.371 1.831 0.704 3.071 7.540
sigma_lambda_H[1] 8.133 4.391 1.576 7.383 18.032
sigma_lambda_H[2] 8.258 4.638 0.880 7.639 18.192
sigma_lambda_H[3] 6.220 4.053 0.950 5.361 15.930
beta_H[1,1] 6.565 1.308 3.432 6.780 8.486
beta_H[2,1] 9.736 0.620 8.373 9.766 10.875
beta_H[3,1] 7.829 0.779 6.006 7.932 9.094
beta_H[4,1] 10.170 7.661 -5.582 10.175 25.535
beta_H[5,1] -0.304 2.438 -5.456 -0.164 3.989
beta_H[6,1] 3.350 4.177 -6.915 4.920 8.307
beta_H[7,1] 1.125 5.542 -10.547 1.417 11.477
beta_H[8,1] 1.215 3.163 -2.633 1.215 3.497
beta_H[9,1] 13.692 5.738 2.770 13.421 25.854
beta_H[10,1] 6.861 1.815 3.132 6.923 10.413
beta_H[11,1] 4.884 3.547 -2.649 5.397 9.953
beta_H[12,1] 2.611 1.103 0.591 2.531 5.114
beta_H[13,1] 9.092 0.858 7.288 9.145 10.521
beta_H[14,1] 2.125 1.041 0.010 2.130 4.185
beta_H[15,1] -6.967 3.451 -13.497 -7.046 0.123
beta_H[16,1] 3.921 3.071 -1.080 3.511 11.277
beta_H[1,2] 7.783 0.266 7.253 7.782 8.311
beta_H[2,2] 9.879 0.166 9.544 9.880 10.206
beta_H[3,2] 8.800 0.204 8.384 8.804 9.197
beta_H[4,2] 3.063 1.441 0.351 3.012 6.003
beta_H[5,2] 1.748 0.985 -0.205 1.766 3.631
beta_H[6,2] 6.103 1.085 3.537 6.294 7.854
beta_H[7,2] 2.422 1.027 0.551 2.357 4.538
beta_H[8,2] 3.032 0.978 1.494 3.160 4.248
beta_H[9,2] 3.217 1.130 1.143 3.193 5.580
beta_H[10,2] 8.025 0.362 7.265 8.038 8.706
beta_H[11,2] 9.827 0.644 8.840 9.730 11.231
beta_H[12,2] 3.937 0.373 3.228 3.926 4.732
beta_H[13,2] 9.146 0.240 8.700 9.137 9.631
beta_H[14,2] 4.017 0.344 3.367 4.011 4.692
beta_H[15,2] 11.529 0.625 10.232 11.533 12.713
beta_H[16,2] 4.448 0.821 2.974 4.410 6.179
beta_H[1,3] 8.437 0.267 7.949 8.423 8.994
beta_H[2,3] 9.858 0.149 9.568 9.858 10.148
beta_H[3,3] 9.460 0.164 9.162 9.452 9.799
beta_H[4,3] -2.183 0.861 -3.851 -2.193 -0.549
beta_H[5,3] 3.900 0.611 2.667 3.904 5.100
beta_H[6,3] 8.334 1.180 6.656 8.012 10.829
beta_H[7,3] -2.572 0.623 -3.796 -2.561 -1.338
beta_H[8,3] 5.268 0.477 4.675 5.207 6.328
beta_H[9,3] -2.745 0.727 -4.196 -2.711 -1.337
beta_H[10,3] 8.550 0.273 8.013 8.552 9.085
beta_H[11,3] 8.504 0.293 7.882 8.533 9.025
beta_H[12,3] 5.268 0.319 4.535 5.304 5.781
beta_H[13,3] 8.866 0.176 8.507 8.871 9.210
beta_H[14,3] 5.789 0.272 5.187 5.806 6.264
beta_H[15,3] 10.290 0.300 9.704 10.288 10.922
beta_H[16,3] 5.902 0.575 4.732 5.917 6.943
beta_H[1,4] 8.096 0.182 7.711 8.101 8.432
beta_H[2,4] 9.851 0.129 9.587 9.853 10.096
beta_H[3,4] 9.951 0.152 9.613 9.963 10.217
beta_H[4,4] 11.821 0.466 10.882 11.822 12.733
beta_H[5,4] 5.574 0.731 4.347 5.496 7.245
beta_H[6,4] 7.229 0.929 5.252 7.467 8.535
beta_H[7,4] 8.551 0.339 7.869 8.557 9.203
beta_H[8,4] 6.727 0.243 6.266 6.738 7.146
beta_H[9,4] 7.434 0.448 6.561 7.428 8.286
beta_H[10,4] 7.526 0.220 7.108 7.520 7.979
beta_H[11,4] 9.515 0.194 9.130 9.514 9.903
beta_H[12,4] 7.255 0.198 6.868 7.256 7.668
beta_H[13,4] 9.062 0.138 8.786 9.066 9.334
beta_H[14,4] 7.834 0.218 7.409 7.833 8.279
beta_H[15,4] 9.462 0.226 9.015 9.461 9.914
beta_H[16,4] 9.433 0.237 8.983 9.432 9.897
beta_H[1,5] 8.897 0.151 8.584 8.903 9.190
beta_H[2,5] 10.676 0.102 10.479 10.671 10.889
beta_H[3,5] 10.759 0.164 10.455 10.757 11.092
beta_H[4,5] 8.328 0.492 7.390 8.319 9.323
beta_H[5,5] 5.206 0.600 3.813 5.289 6.167
beta_H[6,5] 8.965 0.656 7.993 8.810 10.398
beta_H[7,5] 6.573 0.338 5.939 6.565 7.257
beta_H[8,5] 8.246 0.200 7.908 8.233 8.641
beta_H[9,5] 8.189 0.456 7.278 8.185 9.125
beta_H[10,5] 10.040 0.229 9.577 10.038 10.484
beta_H[11,5] 11.585 0.224 11.150 11.586 12.027
beta_H[12,5] 8.608 0.186 8.241 8.603 8.977
beta_H[13,5] 10.133 0.127 9.883 10.130 10.396
beta_H[14,5] 9.242 0.221 8.836 9.230 9.695
beta_H[15,5] 11.322 0.240 10.861 11.320 11.789
beta_H[16,5] 9.855 0.189 9.461 9.859 10.210
beta_H[1,6] 10.119 0.232 9.734 10.088 10.709
beta_H[2,6] 11.472 0.118 11.240 11.472 11.700
beta_H[3,6] 10.620 0.155 10.284 10.632 10.897
beta_H[4,6] 12.940 0.872 11.178 12.943 14.609
beta_H[5,6] 5.662 0.628 4.560 5.622 6.982
beta_H[6,6] 8.754 0.571 7.478 8.797 9.815
beta_H[7,6] 10.556 0.588 9.369 10.572 11.661
beta_H[8,6] 9.534 0.258 9.050 9.543 9.950
beta_H[9,6] 8.291 0.765 6.768 8.295 9.863
beta_H[10,6] 9.140 0.347 8.389 9.162 9.753
beta_H[11,6] 10.703 0.363 9.939 10.729 11.353
beta_H[12,6] 9.384 0.239 8.940 9.382 9.881
beta_H[13,6] 11.092 0.153 10.809 11.088 11.417
beta_H[14,6] 9.916 0.272 9.357 9.916 10.441
beta_H[15,6] 10.523 0.423 9.696 10.525 11.359
beta_H[16,6] 10.576 0.267 9.995 10.587 11.075
beta_H[1,7] 10.202 1.198 7.137 10.465 11.852
beta_H[2,7] 12.113 0.524 11.020 12.145 13.071
beta_H[3,7] 10.289 0.654 8.876 10.340 11.405
beta_H[4,7] 1.315 4.492 -7.488 1.276 10.179
beta_H[5,7] 5.822 1.982 1.905 5.840 9.994
beta_H[6,7] 7.331 3.157 -0.761 8.468 11.122
beta_H[7,7] 7.623 2.991 2.029 7.601 13.805
beta_H[8,7] 10.858 0.992 9.219 10.853 12.386
beta_H[9,7] 4.954 3.884 -3.111 5.128 12.683
beta_H[10,7] 10.062 1.623 7.211 9.940 13.545
beta_H[11,7] 11.345 1.762 8.192 11.175 15.170
beta_H[12,7] 9.755 1.010 7.438 9.842 11.443
beta_H[13,7] 11.721 0.728 10.133 11.798 12.948
beta_H[14,7] 10.649 0.919 8.793 10.678 12.419
beta_H[15,7] 14.060 2.222 9.677 14.035 18.329
beta_H[16,7] 12.927 1.437 10.436 12.813 16.000
beta0_H[1] 8.747 16.656 -26.111 8.984 42.368
beta0_H[2] 10.522 7.986 -6.174 10.664 26.144
beta0_H[3] 9.463 10.515 -11.638 9.584 30.777
beta0_H[4] 3.899 182.336 -365.263 4.762 364.242
beta0_H[5] 4.374 28.665 -51.713 4.267 60.587
beta0_H[6] 5.450 52.572 -119.177 7.449 119.232
beta0_H[7] 3.307 141.312 -276.076 0.944 295.841
beta0_H[8] 6.008 25.978 -16.240 6.407 30.784
beta0_H[9] 7.324 119.152 -231.545 9.138 245.560
beta0_H[10] 8.059 36.370 -64.613 9.039 84.980
beta0_H[11] 8.736 51.426 -103.857 9.240 113.849
beta0_H[12] 6.808 12.276 -18.700 6.683 31.847
beta0_H[13] 9.789 11.732 -12.766 9.980 31.653
beta0_H[14] 7.327 11.086 -15.522 7.077 29.902
beta0_H[15] 6.940 119.782 -244.257 9.296 248.620
beta0_H[16] 7.338 31.493 -56.329 7.462 75.334